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Abstract 



A new approach to formulate the fermion field on lattice is introduced by 
proposing a new Dirac operator on lattice. This approach can eliminate the 
Fermion doubling problem, preserve the chiral symmetry and get the same 
dispersion relation for both Fermion and Boson fields. Then the Weinberg- 
Salam model on lattice may be formulated in this approach. 
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In lattice field theory, it is well-known that the naive discretization of Fermion field 
suffers from the problem of Fermion doubling, while the scalar field doesn't. The studies of 
Fermion species doubling problem persists a long duration since Wilson |]J . To cure from 
this Wilson [|TJ [0 added a new term, so-called the 'Wilson term", to the lattice Fermion 
Lagrangian. This term kills the superfluous components of the Fermion field but breaks 
the chiral symmetry and leads to different dispersion relations for Fermions and Bosons. 
Kogut and Susskind [[| f§] proposed the "stagger model", but it needs four generations of 
Fermions at least, and deals with the Fermions and Bosons on different footing. Besides 
these two popular prescriptions, Drell et al || developed another approach which preserves 
the chiral symmetry and at the same time correctly counts the number of Fermion states. 
Their crucial point was to introduce a lattice gradient operator which couples all lattice 
sites along a given direction instead of coupling only nearest-neighbor sites. However, no 
convincing approach which can fit all the physical phenomena exists so far |J, such as we 
don't know how to reformulate the Weinberg-Salam model on lattice. 

In this paper, to overcome those difficulties we propose to adopt a new approach in which 
a new Dirac operator is introduced. The square of this Dirac operator is the D'Alembertian 
operator, which recovers the relation in continue field theory. By using this approach, we can 
eliminate the Fermion doubling problem and at the same time preserve the chiral symmetry. 
Besides these, we can get the same dispersion relation for both Fermions and Bosons. 

First let us take a look upon how the Fermion-doubling problem arises. For free scalar 
field (fi(x), the lattice action in the Euclidean space is (for more detail, see ) 



(1) 



With the help of displacement operator which is defined as 
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the free scalar action ([[]) may be written as 
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where the discrete derivative = — 1 and D'Alembertian operator □ = J2n 
Using the formula 
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we can transfer (ffl) into the momentum space form 
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(2 

The dispersion relation read from this equation is 



W = E^in 2 (^) + M 2 . (6) 



The naive dicretization of free fermion field is usually taken as 
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S(1>) = Ely E - Wn + Ma 4 ^i} 

n /•( 
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= E{ V E - ^-/^n + Ma 4 ii} (7) 

where Tp = 7^(9^ is the Dirac operator. Transferring it into momentum space we get 

SW = I ^H)[^7^ + M]^W. (8) 
J (2tt) 4 ^ a 

From this equation we see the dispersion relation 

S4k) = Y / ^^ 2 (ak fl )+M 2 (9) 
n a 

is different from (|6]), and the doubling of Fermionic degrees of freedom appears. 

In above formulas, two key operators on lattice were introduced. The first is the Dirac 
operator 'p, the second is the D'Alembertian operator □, which domain the dynamics of 
fields. Now we may raise a question, why the lattice description of scalar field is successful, 
while the fermion field fails? One possible way to answer this question is to regard the 
lattice definition of D'Alembertian operator is "good" but that of Dirac operator is "bad". 



So we advocated a new point that all the disaster of lattice Fermions come from the "bad" 
definition of Dirac operator. Therefore, to overcome the Fermionic difficulties we propose a 
new lattice Dirac operator, which may shed the new light on studying lattice Fermion fields. 
Consulting the definition of the translation operator in the momentum spaceQ. 

dk 4 
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a Jkna iP(k)e ik » a (10) 
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we can define an operator Re , corresponding to the "half-spacing translation" , as 



R^ n := ^ J^y^me 
R-^n ■= Jl ^e ik ™me-^- (11) 



From the definition ( [TlD we can also see that Re and R_e are not "local operators" 



When acting on some functions of lattice sites, they not only concern with the nearest 
neighbor sites but also all the sites in direction fi on the lattice. However, it is evident that 

R /j- —R_ m 

the combination 7 - — 21 will approximate to when the lattice spacing a — ► 0. So we 

R j-i — R_ P" 

can call the combination 7 - — — the " quasi-local operator" . The quasi-local feature is the 
price for getting the above mentioned nice characters. 

In this approach, the key point is to introduce the auxiliary operator Re. Actually, 
physical fields (p and ip distribute on the integer sites and the elementary translation operator 
is R^. To discuss the physical meaning of what we really did, it is desirable to represent the 
auxiliary operator Re in terms of physical operator R^. 

By the help of definition (|TT|), it is easy to show 



Z7T m=-oo J — 

-\ oo I -r \n— m 

E 1 '. R^n. (12) 



m=— oo 2 



1 The similar operator was also introduced in some papers at different angle H| 
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So in terms of R^, the auxiliary operator Re may be written as 



2 

OC 



Re = - V ^4=- i£» . (13) 

" m=-oo " t ^ 2 



Similarly, 



R-^ = ~- E ^T^V (14) 

From the above formulas, we see that the operators.^ and R_e are well defined globally, 
although they were introduced in the first Brillouin zone. 

Using the expressions of operators Re, R-e in equations (|13|) and ([[4]) we may prove 
the following relations hold, 

Re = R„, (15) 
ReR—e — R—eRe — I- (16) 

2 2 2 2 v ' 

If the Dirac operator is redefined as 

Re — Re 

P = 7m( 2 fl 2 ): (17) 
by the help of equation fll3] ) and (|14|) , it may be written as 

? = — E L 4 L T E i»{R-» - K)- (i 8 ) 

It is interesting to find that in this case we may prove the following equation exits 

U=ip\ (19) 

which recover the same relation of continuum field theory. In continuous limit a — > 0, 
□ — > Ylt=i 9^, because the square root of D'Alembertian operator is unique in continuum 
theoryso by the help of equation (|TJ|), we show that the continuous limit of the new Dirac 
operator is just the ordinary Dirac operator. 

Substitute the modified Dirac operator into (|7]), we may get the free fermion action as 

9 00 ( — -\\ m m 

Sty) = E E &7„(#?„ - + Ma 4 ^ n }. (20) 

n 71 771=1 m 4 n 



which shows how the fermion field on non-neighbour sites correlated each other. When m 
is taken to 1, (|20 ) reduce to the form of naive discretization action ([7p. Transform into 



its momentum space form, we have 

d A k , vr ^ , 4z N ~ (-l) m msm(mak 



/" d Ac , vrV -^ , 4z. ^ -1 msin mafcJ „ , .,. ,„ . 

- / j^H-k) E 7 ,(~) E J (m ,_l )a + "hW (2D 

= / ^(- fc )^E7^ + M]^(fc), (22) 



(27T) 

where function A fx) is defined as 

sin (| — pir), (2p — lW < x < (2p + lW 

A(x) = { K2 V . (23) 

0, x = (2p+l)n 

and may be illustrated by following figure, 
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From this we see immediately that there is no double- counting anymore and the Fermion 
field has the same dispersion relation as the one for scalar field (^j) one on the boundaries of 
the Brillouin zone A(a/c M ) = 0. When M = this new formula ( pop will preserve the chiral 
symmetry. 

The next work is to construct the gauge invariant action in this frame. Under gauge 
transformation of lattice theory, will be changed into its covariant form as 

R^n > U n ^R^ n 

R-^ n ► Ul-wR^n (24) 

where U n ^ are the known link variables. Then we have 

R™^n — {Ul-^R-^n (25) 

Expanding the exponential, one gets 



(Un,fj,Rfi) Un^Un+n^ • • • fji+(m-l)]tt,/i-^u 

(jA C> \m _ tj\ tj] . . . jj Dm /of?\ 
\ w n—/j,,fj, ± ''—nJ u n—fi,fj,'- y n—2[i,fi w n— m/i,^- 1 L — v z ' u / 

Now we can build a gauge invariant toy model including a Fermion field if), scalar field 

(f) and gauge field U n ^, where U njfl is the link variable of gauge group G. The action of the 
toy model may be written out as following 

S Sp erm i on -\- Sgcalar Sg aU g e (27) 



where 



2 00 (—'\) m rn 

S Fermi on = 5> S - E W^H^nl^Ul^r - (WM 
n 71 m =l m 4 (i 



2 4 

Sscalar = YA^^nJ^ + U^R^ - 2]0 n ] 2 + ^ 2 ^}, (28) 



" Plaqu 



where tjj L and ip R are the left and right hand Fermion respectively. 

With all the previous points in mind, a construction for the SU(2) x U(l) electroweak 
theory on lattice is now available. Consider just the leptonic first generation section: electron 
and electron neutrino and define the fields as follows: 

M 

L — , R = e R 

\e L ) 

U r l, = U n ,,V„U^ = V^ (29) 

where U n>fl exp l9aT ' An ^ is SU(2) gauge field and V n ^ = exp~2 l9 ' aBn ^,V' = exp"' 9 '"- 8 "'" are 
U(l) gauge field. Because the action of scalar field and gauge field parts are easy to write 
out, we only give the action of fermionic part here, 

2 00 (— l) m m - L L 

SFermion = E{° ~ E 2 _ I E [(C^-R-^)™ ~ (U n ^R^) m )L n 

n 71 m =l rn 4 n 

+E{« 3 - E yrEWW" - (U^R,r)R n 

n n m =l 171 4 H 

+\a\L n <p n R n + R n <PiL n )}, (30) 



where A is Yukawa coupling constant and 



is the Higgs doublet. 



In summary, in this paper we have proposed to introduce the Dirac operator Tp = 

00 ( — \\m m 

^ E — 2 r(-^-/i — R™) m lattice theory, which shares the property^* 2 = □. By doing 

m =\ m ~ 4 

so, we can eliminate the Fermion-doubling problem, get the same dispersion relation for 
both Boson and Fermion fields and preserve the chiral symmetry. We can also construct 
the gauge theory in the new frame. Because of these good characters, we can discuss chiral 
model which is difficult to touch upon in normal lattice theory and we completed the con- 
struction of Weinber-Salam Electro- Weak model on lattice. Modifying the Dirac operator is 
not new in lattice theory, actually, the introduction of Wilson fermion may also be regarded 
as another kind of modifying method in some sense. 
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